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■ Abstract: Contributions of liadronic effects to tlie muonium physics and anomalous magnetic 

' moment of muon are considered. Special attention is paid to higher-order effects and the 

fT^ ' uncertainty related to the hadronic contribution to the hyperfine structure interval in the 

ground state of muonium. 
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[Traduit par la redaction] 



1. Introduction 



o 

■ Quantum electrodynamics (QED) is a theory which covers all interactions of leptons (electrons and 

I muons) with photons. However, it may be not sufficient even for pure leptonic systems that are not free 

of hadronic effects appearing because of hadronic intermediate states. Such effects cannot be calculated 
il ab initio and additional data on these states are needed. Since the data can be achieved mainly from 

1^ ■ experiment, their availability and uncertainty impose some principal limits on any QED calculations. 

Recently a number of projects on intensive muon sources have appeared [ 1 ] and a new generation of 
precision experiments on muonium may appear in near future. Therefore it is timely to investigate the 
^> ] principal limits of QED tests related to hadronic physics. Here we discuss hadronic effects in muonium 

5_i ■ and free muon. 

2. Hadronic effects in muonium hyperfine structure: leading contribution 

The leading hadronic contribution to the hyperfine splitting of the muonium ground state comes 
from the diagram depicted in Fig. 1 . We separate QED and hadronic effects in this diagram [2] 

2 E f 

Ai^hadr (leading) = -2— " ^ / dsKMu{s)p{s) (1) 

TT-^ mf^ - mi 1 + ttfj, J 
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Fig. 1. The leading hadronic contribution to the muonium hjfperfine structure 




Fig. 2. The skeleton two-phonon exchange diagrams 




and evaluate them independently. Here Kmu{s) is a result of the calculation of a two-photon exchange 
diagram with one massive photon (s = A^), the Fermi energy 



E 



8 4 mlml 
F = -{ZaY— — , ^ ^3 (1 + a^) 



(2) 



is a nonrelativistic value for the hyperfine structure interval, s is the center-of-mass energy squared and 
p{s) is the hadronic spectral function: 



2 , , cr(e+e — > 7 ^ hadrons) 



(3) 



The Fermi energy determined above contains the anomalous magnetic moment of muon, a^, which is 
also affected by hadronic effects (see. Sect. 4). 

The basic expression has been obtained in Ref. [2]. Here we briefly reproduce the derivation. We 
start from the evaluation of the two-photon skeleton contribution (Fig. 2) 



AE'skel = Ep 

where 
J(fc,M) = 2 



6 a TO^TTie 

TT ml -ml 



J -^[J{k,m^) - J{k,me)j 



(4) 
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We introduce hadronic polarization effects with a dispersion substitute for the photon propagator 
1 a f ds p{s) 



fc2 TT ./ fc2 -j- S ' 

and after fc-integration we obtain 



(5) 



^Mu(s) = 



4ml 



1 



4ml 1 
!^ln- 



4m„ 



4mf, 



S _|_ S 

4ml 
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(6) 
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We are interested in relatively high s and the expansion of the exact result 



Kmu{s) 



Ami 



9 , s 15 



^ In * '^'^ 



16 



96 



+ 



Ami 



21 , s 

In — 

128 ml 



73 
512 



+ ... 



(7) 



is useful for further evaluations. For s 



the leading term 



15 
16 



(8) 



is a good enough approximation (within 2%). For higher .s it is even better (see Fig. 3). 

A simple model can be used to estimate the leading hadronic contribution (cf. Ref. [3]). The spectral 
function consists of several narrow-pole contributions (p, w, (j)) 



Ppole(s) 



" res 



res) 



37rr(res ■ 



Mr, 



5{s - M^J 



(9) 



and background with i? = 2 for 1 Gev^ < s < 4 Gev^ and = 4 for 4 Gev^ < s. Here Mrcs, r(res 
e+e") and /res are the resonance mass, leptonic width and coupling constant, respectively. The results 
are summarized in Table 1 . We estimate a systematic uncertainty of the model as 1 0% . That is not good 
enough for a calculation of the leading term but sufficient for higher-order hadronic corrections. 
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Table 1. Estimation of the leading hadronic contribution to the hyperiine structure in a simpUfied model. 



Contribution 


AS 






[kHz] 


[%] 


P 


0.151(7) 


65(3) 




0.013 


5.5(2) 


<A 


0.014 


6.0(2) 


s = 1 - 4 Gev^ 


0.045 


19 


s > 4 Gev^ 


0.009 


3.8 


Total 


0.232(7) 


99(3) 



Any realistic calculation of the leading hadronic term targeting a higher accuracy should deal with 
experimental data on e+e^ annihilation into hadrons (see e.g. Refs. [4, 5, 6]). The most accurate recent 
result [5]' is 

Ai;hadr (leading) = 0.233(3) kHz . (10) 

Since the publication of Ref. [5], some new experimental data appeared and the result can be now 
updated. These new data are: 

• The CMD-2 detector has recently published the final results of their analysis for the reaction 

e+e" TT+TT" in the c. m. energy range 610 to 960 MeV [7]. Because of the new, more precise 
approach to calculating radiative corrections which also include a correction for the final state 
radiation, the values of the cross section are slightly smaller than before. 

• The new data on the (f) meson leptonic width obtained by the CMD-2 and SND groups in Novosi- 
birsk [8]. 

As a result, the updated value of the leading order hadronic contribution becomes 231.2 ± 2.9 Hz 

compared to our last year value of 233.3 ± 3.1 Hz quoted above. We expect that after the analysis of 
other hadronic modes is completed in Novosibirsk, the precision of the leading order contribution can 
be improved to about 2 Hz. 

3. Hadronic effects in muonium hyperfine structure: next-to-leading term 

With an accuracy at a one-percent level for the leading hadronic contribution one has to take into 
account higher-order hadronic effects. The corresponding diagrams for some higher-order hadronic 
terms are presented in Fig. 4. The first evaluation presented in Ref. [2] gave a result 

A£hadr(non - leading) = 0.007(2) kHz . (11) 

Most of the contributions can be expressed in terms of corrections to i^Mu in the form 

a ^ml (9 ml s 9 o s 15 mf. s 
AK = Jjln^-|ln— + ^ln2 — + -|ln-|+Ciln— + C2 

a 4:1711 

~ ^{47.6-|-17.7-M0.0-|-(7i ■4.O-I-C2} . (12) 

TT S 

^ The result of Ref. [6] has been obtained using r lepton data and the accuracy quoted there is overestimated. The discussion 
of the possibility to use r data can be found in Sect. 4. 



???? NRC Canada 



Eidelman et al. 

Fig. 4. The next-to-leading hadronic contribution to the hyperfine splitting in muonium 







Table 2. Higher-order hadronic contributions to muonium h3fperfine splitting. The results are presented in units 

of a/7r • Bhadr (leading). 



Correction 


Contribution 




[a / - ■ -t'hadr] 




9.8 




1.6 (7) 


ASc 


-1.0(3) 




±0.1 


ASe 


-1.0(3) 


AEt 


±0.1 


AEg 


0.10(2) 


AEh 


±0.1 


AEi 


±0.3 


ABhigh 


9.2(14) 



The numerical values are given for s = m^. We present here the calculation of double-logarithmic 



terms and that of a single-logarithmic term with the bigger logarithm ln(m^/mg'' 
4. A calculation of other logarithmic contributions is in progress. 
The biggest contributions come from Fig. 4a. The result is 



11 > Iniml/ml) 



3 a 



mf, 



3w s 
Our final estimate is 
A£;hadr(non - leading) = 0.005(2) kHz 



27 
32 



+ 



■In 



s 15 
m?. 16 



2 

mf. 
In -I 



(13) 



(14) 



The contributions are summarized in Table 2. The result is somewhat lower than in Ref. [2] mainly 
because of single-logarithmic contributions from Fig. Ac and e. 

Recently higher-order hadronic effects were also studied in Ref. [9]. In particular, the hadronic 
light-by-Ught contributions were under examination. The contribution of the diagram in Fig. Ah is 
1.5 mHz and so it is consistent with our conservative estimation. A bigger contribution of 3.65 Hz 
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Fig. 5. Leading hadronic contribution to the anomalous magnetic moment of muon 




found there (related to Fig. 2 of Ref . [9]) is actually due to the correction to the anomalous magnetic 
moment of muon and affects hyperfine spUtting accordingly to Eq. (2). Therefore it should be excluded 
from our consideration. 

More accurate calculations are in progress and we hope to reduce the uncertainty by a factor of 4 
[10]. 



4. Anomalous magnetic moment of muon 

Since the 2001 pubUcation of the E821 result on the anomalous magnetic moment of muon [11], 
a lot of efforts have been put on the new estimation of the hadronic correction to the muon anomaly, 
ghad^ see e.g. Ref. [12]. However, these works were based on either not final experimental data or 
involved various theoretical assumptions. In view of the extreme importance of the E821 result and 
possible indications to deviation from the Standard Model it is necessary to have a model-independent 
analysis based exclusively on the data, similar to that of Ref. [13]. The new calculation based on the 
most recent data from e+e" annihilation at low energy in Novosibirsk and Beijing as well as from 
r decay data, which also includes new developments in theory, has recently been pubUshed [14]. It 
showed a real breakthrough in the accuracy of the a^"'' estimation which became possible after the 
final data on e+e" ■k^it~ from CMD-2 with a 0.6% systematic uncertainty have been published [7] 
and, in addition, the final analysis of the R measurement by BES between 2 and 5 GeV appeared [15]. 
The resulting accuracy of the estimate increases dramatically and reaches ~ 70 ■ 10"^'^ based on e+e— 
data only. This accuracy is much better compared to the uncertainty of 153 • 10~^^ obtained in 1995 in 
Ref. [13] and is very close to that obtained in the 1998 analysis of Ref. [16] which also included the r 
data and QCD sum rules. Calculations show that another significant improvement (by a factor of about 
1.5) is possible when high precision r lepton decay data are also used. However, for their reUable use 
one should first understand the reason of the observed discrepancy between two data sets: the two pion 
spectral functions obtained from e+e~ and r lepton data do not agree confirming earlier evidence for 
this effect [17]. The solution of this problem should also involve a thorough investigation of the effects 
of isospin breaking corrections as well as additional radiative effects in r decays [18]. 

The real uncertainty of the theoretical estimate will also strongly depend on the understanding of 
the role of higher-order effects. These include diagrams with two loops which can be reliably calculated 
fromthee+e~ data using dispersion relations [19] and much more complicated effects due to hadronic 
light-by-Ught scattering. At the present high level of accuracy, both effects give a quite substantial 
contribution of the order of 100 • 10~^^ each, i.e. about or even larger than the current uncertainty 
of the leading order hadronic contribution. Recent reevaluations of the pseudoscalar pole contribution 
dominating in the hadronic light-by-light scattering diagrams revealed that its sign was previously 
wrong [20]. This initiates new calculations of the higher-order effects and more realistic estimates of 
their accuracy [21]. 
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5. Summary 

The hadronic effects become an unavoidable part of theory of "pure" QED quantities, such as 
the anomalous magnetic moment of muon and hyperfine interval in the ground state of muonium. In 
the case of muon g-2 they limit present accuracy of the theoretical prediction whereas for muonium 
their impact is more relevant for the future experiments with intensive muon sources and thus higher 
expected accuracy. 

The new experimental result on the anomalous magnetic moment of muon has recently appeared 
[22]. It has about the same accuracy as theory values from e+e" and/or r data. A more accurate result 
is expected soon. The hadronic contributions dominate in the theoretical uncertainty and require a 
more detailed study. On the contrary, no results on muonium hyperfine structure are expected in near 
future and theory is more accurate than the experiment. However, to clearly understand what level of 
accuracy could be a target of new experiments on muonium, one should study theoretical uncertainties. 
The bound state QED theory is overviewed in Ref. [5] and some progress for ab initio calculations is 
possible. Better understanding of hadronic effects is needed to find a level of accuracy which cannot be 
superseded by the ah initio QED calculation. Several projects on the calculation of hadronic effects in 
physics of muon and muonium which are shortly described here are in progress and we hope to report 
on their results soon. 
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